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The 5.4 kton MINOS far detector has been taking charge-separated cosmic ray muon data since the 
beginning of August, 2003 at a depth of 2070 meters-water-equivalent in the Soudan Underground 
Laboratory, Minnesota, USA. The data with both forward and reversed magnetic field running 
configurations were combined to minimize systematic errors in the determination of the underground 
muon charge ratio. When averaged, two independent analyses find the charge ratio underground to 
be 

^M+M.- = 1-374 ± 0.004 (stat.)t°:;]lg (sys.). 
Using the map of the Soudan rock overburden, the muon momenta as measured underground were 
projected to the corresponding values at the surface in the energy range 1-7 TeV. Within this 
range of energies at the surface, the MINOS data are consistent with the charge ratio being energy 
independent at the two standard deviation level. When the MINOS results are compared with 
measurements at lower energies, a clear rise in the charge ratio in the energy range 0.3 - 1.0 TeV is 
apparent. A qualitative model shows that the rise is consistent with an increasing contribution of 
kaon decays to the muon charge ratio. 



PACS numbers: 

I. INTRODUCTION 

The MINOS far detector is a 5.4 kton calorimeter, with 
magnetized steel planes located at the Soudan Under- 
ground Laboratory, Minnesota, USA. It was designed to 
study neutrino oscillations with the NuMI beam which 
originates 735 km away at Fermilab. At a depth of 710 
meters below the Earth's surface, the MINOS far detector 
is the deepest experiment to measure cosmic ray muons 
with a magnetized detector, thus providing a capability 
to distinguish /i"*" from /i^ with large statistics. The data 
correspond to muon energies at the surface in excess of 
approximately 1 TeV. Above this energy, there are few 
measurements of the charge ratio /N^- . 

Cosmic ray muons are produced when primary cosmic 
ray nuclei, mostly single protons, interact at the top of 
the atmosphere to produce hadronic showers. The pions 
and kaons in these showers decay to muons, which are 
measured in detectors on the surface and underground, 
and with balloon-borne experiments in the atmosphere. 
Since the cosmic ray primaries are positively charged, 
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there are more positive than negative pions and kaons 
in the resulting hadronic showers. In the energy range 
from 3 to 100 GeV, the CORT (Cosmic- Origin Radia- 
tion Transport) cosmic ray Monte Carlo [l[ predicts a 
relatively constant charge ratio of N^+/N^- ss 1.3. In 
a compilation of measurements, the charge ratio was ap- 
proximately constant at N^+/N^- k, 1.27 with uncer- 
tainties increasing from ~ 1% at a few hundred MeV to 
^ 6% at 300 GeV [;2|. A more recent measurement by 
the L3-I-C experiment at CERN found similar results in 
the range 20 - 500 GeV [1]. 

At energies greater than a few hundred GeV, several 
competing processes can affect the charge ratio. The 
muons seen by MINOS result from pions and kaons that 
decay before they interact in the atmosphere. As energy 
increases, the fraction of muons seen from kaon decays 
also increases because the longer-lived pions have be- 
come more likely to interact before decaying than the 
shorter-lived kaons. Consequently, kaon decays begin 
to make an increasingly more important contribution to 
the muon charge ratio at these energies. Since strong 
interaction production channels lead to a muon charge 
ratio from kaon decays that is greater than that from 
pion decays, the measured charge ratio is expected to 
increase. Several competing processes, however, could 
counter this increase at even higher energies. Decays of 
charmed hadrons are one such process. There is also 
the possibility that heavier elements become a more im- 
portant component of cosmic ray primaries as the energy 
increases. This increasingly heavy composition would de- 
crease the ratio of primary protons to neutrons, thereby 
decreasing the muon charge ratio. With careful measure- 
ments of the /N^- ratio in the cosmic rays, models 
of the interactions of cosmic rays in the atmosphere can 
be improved. 

In addition, measurements of the cosmic ray muon 
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charge ratio from a few GeV to a few TeV are impor- 
tant to constrain calculations of the atmospheric neutrino 
fluxes. These are of interest both for detailed measure- 
ments of neutrino oscillations in atmospheric neutrino 
experiments and also for calculations of backgrounds for 
neutrino telescopes. The muon charge ratio is a partic- 
ularly useful tool for testing the predicted atmospheric 
vjv ratio 

In MINOS, underground charge-separated cosmic 
muons were first studied in detail by Rebel [BJ and Beall 
0. 



II. THE MINOS FAR DETECTOR 

The MINOS far detector is a steel-scintillator sam- 
pling and tracking calorimeter located at a depth of 
2070 meters- water-equivalent (m.w.e.) in the Soudan 
Underground Laboratory in an iron mine in north- 
ern Minnesota (47.82027° N latitude, and 92.24141° W 
longitude) j3] . The detector is made of 486 octagonal 
planes of 2.54 cm thick steel laminates, interleaved with 
484 planes of 1 cm thick extruded polystyrene scintillator 
strips at a 5.94 cm pitch. Each scintillator plane has 192 
strips of width 4.1 cm. The length of each strip depends 
on its position in the plane and varies between 3.4—8.0 m. 
The scintillator strips in alternating detector planes are 
oriented at ±45° to the vertical. The modular detector 
consists of two supermodules (SM) separated by a gap of 
1.1 m. Fig. [T] shows the coordinate system used in the 
MINOS far detector cosmic ray analysis. In terms of this 
coordinate system, events are described as coming from 
the zenith angle d (the polar angle measured from the y- 
axis) from vertical (0 = 0°) toward the horizon {d = 90°), 
and the azimuthal angle measured in the x — z plane 
from true north (0 = 0°) to the east ((/) = 90°). 

Light from charged particles traversing the MINOS 
plastic scintillator is collected with wavelength shifting 
(WLS) plastic fibers embedded within the scintillator 
strips. The WLS fibers are coupled to clear optical fibers 
at both ends of a strip and are read out using 16-pixel 
multi- anode phototmultiplier tubes (PMTs). The sig- 
nals from eight strips, each one of which is separated by 
approximately 1 m within the same plane, are optically 
summed and read out by a single PMT pixel. The fibers 
summed on each pixel are different for the two sides of the 
detector, which enables the resulting eightfold ambiguity 
to be resolved for single particle events. For all other 
types of events, ambiguities are effectively resolved using 
additional information from timing and event topology. 

The data acquisition and trigger have been described 
in 0. The primary trigger requires activity to be ob- 
served on 4 planes out of 5 within 156 ns. The detector 
calibration has been described in |1Q[. More detailed de- 
tector information can be found in [8|. 

In order to measure the momentum of muons travers- 
ing the detector, the steel has been magnetized into a 
toroidal field configuration. A finite element analysis cal- 



culation of the magnetic field strength for a typical MI- 
NOS detector plane is shown in Fig.O These calculations 
show that each SM is magnetized to an average value of 
1.3 T by the 15 kA current loop that runs through the coil 
hole (c.f.. Fig. [T]) along the detector's z-axis. The field 
is saturated near the coil hole at a strength of approxi- 
mately 1.8 T, falling to almost 1 T near the edges. There 
are small variations in field strength near the corners and 
along the gaps between the eight plates that make up 
each steel plane. The detector field was designed to bend 
or 'focus' negatively charged muons travelling from de- 
tector South (i.e., /i~ resulting from i/^ interactions in 
the detector from neutrinos originating in the Fermilab 
NuMI beam) toward the center of the detector. 

The toroidal magnetic field of the MINOS far detector 
impacts the acceptance of /i^ and /i"*" entering the detec- 
tor as a function of their incoming trajectories and the 
field direction. A muon is "focused" when the magnetic 
field steers it toward the center of the detector and "defo- 
cused" when directed away from the center. These effects 
are most apparent for muons with trajectories that are 
parallel to the detector z axis. In one field orientation, 
Forward Field running, [f which enter the detector from 
the south and the which enter from the north will be 
focused into the center of the detector, while the muons 
in the opposite charge sign and trajectory combinations 
with be defocused. Forward Field running ("DF" - data 
forward) is the default configuration for MINOS data- 
taking with the NuMI beam from Fermilab. MINOS has 
a second field orientation. Reverse Field running ( "DR" - 
data reverse), in which the coil current is reversed and /i"'" 
from the south are focused into the detector. These fo- 
cusing/defocusing effects are most important at the edges 
of the detector acceptance and as a result, the charge ra- 
tio for muons with incoming trajectories on the edges of 
the detector acceptance will be either enhanced or sup- 
pressed depending on the charge and incoming direction 
of the muons. 



III. DATA ANALYSIS 
A. Cosmic Muon Data Sample 

In this paper we present results from data recorded be- 
tween August 1, 2003 and February 28, 2006. During this 
period, the detector ran with both the DF and DR mag- 
netic field configurations. There were 609.82 live days of 
Forward Field running and 201.75 live days of Reverse 
Field running. There were 28,994,380 events in the DF 
sample and 8, 898, 551 events in the DR sample. 

B. Event Selection 

The first stage of the event selection is to identify and 
remove periods of data associated with detector hardware 
problems. The criteria defining bad runs are described 
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FIG. 1: The coordinate system for the MINOS cosmic ray 
analysis. The octagonal steel and scintillator planes of the 
MINOS far detector are 8 m across. The detector is 30 m 
long. The central hole is for the magnet coil. The +z-axis is 
along the long axis of the detector and points toward detector 
North (N). Detector South (S) points back towards Fermilab. 
The j/-axis is directed toward the zenith. The a;-axis direction 
is chosen to make a right-handed coordinate system. The ori- 
gin of the coordinate system is the center of the South face of 
the detector. Detector North (N) is rotated from true North 
by an angle a — 26.5548° about the y-axis as measured by a 
gyro-theodolite; detector North therefore points along an az- 
imuthal angle of 333.4452°. Alternating planes of scintillator 
strips are oriented along either the u or v axis directions, a 
coordinate system in which the x~y plane is rotated by -1-45° 
about the z-axis. 

in The sample was selected using a series of cuts 

that are described in detail below and the numbers of 
events remaining at each stage in the selection are listed 
in Table [J 

1. Pre-Analysis Cuts 

The first cut in the event selection requires at least 
one reconstructed track in the event ("1. no reconstruc- 
tion"). This requirement predominantly removes noise 
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FIG. 2: Finite element analysis model for the toroidal mag- 
netic field in a plane of MINOS steel. The coordinate system 
for the field map is shown in Fig.[T] In this map, the detector 
plane is being viewed from detector North. 



where the primary trigger was satisfied, but there was 
not enough activity to resolve the eightfold ambiguity 
from the optical summing of the scintillator strips. The 
second requires that there is only a single track found 
by the track- fitting algorithm ("2. multiples cut"). The 
third requires that the coil be on and in a known state 
("3. coil status cut"). 



2. Analysis Cuts 

The next set of cuts are meant to separate muon tracks 
from the background with high reliability. These cuts re- 
quire that: a track must cross at least 20 planes in the 
detector ("1. 20 plane cut"); a track must have a path 
length of at least 2 m ("2. 2.0 m track length cut"); the 
entrance point of a track was required to be less than 
50 cm from an outside surface of the detector ("3. Fidu- 
cial cut" ) ; and a track must pass a quality cut based on 
X)itlndf < 1.5 (4. "fit quality cut"). The xju/ndf pa- 
rameter is returned by the Kalman filter [l^l , which is the 
track fitting algorithm used in this analysis [1] . The dis- 
tribution used to select this cut value is shown in Fig. [3] 
and this cut assures that the track found is a good fit to 
the track hit points. 

In Fig. [3] we show the muon rate (Hz) as a function of 
day number from the beginning of data taking with the 
complete and magnetized detector after the pre-analysis 
cuts and analysis cuts. The fluctuations are consistent 
with seasonal variations in the cosmic ray muon flux [l3| . 
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FIG. 3: Distribution of Xfu/'n'd-f for the fits to cosmic muon 
tracks. The cut was made at y^julndj > i.5. The shaded 
region shows the excluded tracks. 
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FIG. 4: Muon rates per day as a function of day number from 
the start of data taking with the complete and magnetized 
detector after the pre-analysis cuts and analysis cuts 1-4. The 
shaded area shows the period of reverse field running. 



Consistent results for the cosmic ray muon charge-sign 
ratio were obtained with the two different approaches. 

(1) "track quahty cut" - the significance of the mea- 
sured muon charge sign and momentum for a track 
was required to be {q/p)/<Jq/p > 2.2, where q is 
the charge sign and p is the fit momentum. Here 
{q/p) is the fit parameter returned by the Kalman 
filter and the uncertainty on this fit parameter is 
<Jq/p- For this cut we treat the quantity {qlp)loq/p 
as positive definite. In Fig. [5] we show A^^+ /N^- , 
the ratio of the number of fj,'^ to as a func- 
tion of {q/p)/crq/p after all other cuts have been 
made since the final two cuts in the analysis are 
closely correlated. For {q/p)/crq/p > 2.2, the charge 
ratio becomes asymptotically flat, suggesting that 
the charge sign and momentum are well fit. As 
{q/p)/(Jq/p tends to zero, the fitter is becoming less 
reliable at determining the charge sign. For the 
lowest values of the significance, the fitter picks 
the two charge signs with equal probability and the 
measured charge sign ratio tends to unity. Events 
with low values of the charge sign significance are 
typically high momentum tracks (> 100 GeV/c) 
that traverse the detector in such a way as to bend 
only slightly in the magnetic field. These events 
have random charge sign identification. 

Fig. m shows that the DF and DR data sets appear 
to behave very similarly with respect to the entire 
suite of charge sign quality cuts. Not withstanding 
this similar behavior, there is a clear offset in the 
charge ratio on the plateau regions where the cuts 
define the tracks with well-determined charge sign. 
For the DF data set, iV^+ /N^,- = 1.40; for the DR 
data set, /N^^- = 1.33. This overall systematic 
difference is discussed below. 



3. Charge Sign Quality Cuts 

For the analysis presented in this paper it is neces- 
sary to cleanly identify the charge sign of selected muons 
and to ensure that systematic uncertainties in this iden- 
tification are minimized. This clean selection is achieved 
by placing charge-sign quality requirements on the recon- 
structed tracks. There are two components to the charge- 
sign quality cuts. The first cut assures that the charge 
sign and momentum returned by the track-fitting algo- 
rithm has been well-determined. MINOS uses a Kalman 
filter technique for track-fitting fl^l that simultaneously 
determines which hits belong on a track and the momen- 
tum of the particle. The technique involves a series of 
recursive matrix manipulations to specify the trajectory 
of the particle as well as the ratio of its charge to its 
momentum. The second cut minimizes residual system- 
atic uncertainties that remain in the data set after the 
cuts already described have been applied. To increase 
the robustness of the result two different approaches to 
identify a clean sample of events were adopted (2a, 2b). 



To check the consistency of the analysis, the final cut 
used was either (2a) or (2b). 

(2a) "MIC" (Minimum Information Cut) analysis - a 
track was required to have at least 60 planes where 
the hit information was within 3.5 meters of the 
detector center; it was also required that the track 
fitter use the hits from these planes in its determi- 
nation of {q/p)- Here a 'plane of hit information' 
is defined as a plane containing a strip from which 
signal is read out on both ends. This cut was moti- 
vated by the need for high quality track information 
to resolve ambiguities in the charge sign determi- 
nation for events with p ^ 50 GeV/c. 

The effect of this cut is shown in Fig. [HI where the 
muon charge ratio, /N^- is plotted for the (a) 
DF and (b) DR data sets as a function of fit mo- 
mentum. In this plot the number of planes of infor- 
mation required where the hits are within a radius 
of 3.5 m from the detector center has been varied 
from (no cut) -100 planes. Without the MIC cut. 
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FIG. 5: The A'^^+ /N^~ ratio for reconstructed muon tracks as 
a function of (q/p) /uq/p after all other cuts, (a) Data Forward 
(DF) data distribution; (b) Data Reverse (DR) data distribu- 
tion. As indicated, the cut was chosen at {q/p)/aq/p > 2.2 for 
both data sets, where {q/p)/o'g/p becomes asymptotically flat. 
The figures are labeled with the flts to a constant 7V^+ /N^- 
ratio for the cut at {q/p)/aq/p > 2.2. The fits are superposed 
onto the data as horizontal lines. 
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FIG. 6: The muon charge ratio Af,j+ /A*'^- for (a) the DF and 
(b) the DR data sets as a function of fit momentum as the 
number of planes of track information within 3.5 m of the 
detector center is varied from zero (no cut) - 100 planes. 



there is a large bump in the charge ratio distri- 
bution that peaks in the neighborhood of 40 — 50 
GeV/c and which decreases as the fraction of the 
track information in the inner part of the detector, 
where the field strength is well-characterized, in- 
creases. A second feature, which has much poorer 
statistics, appears near 100 GeV. If these features 
were real physical effects, then we would not expect 
them to disappear as the quality of the fit informa- 
tion improved nor would we expect them to become 
reversed when the field was reversed. But as is clear 
from Fig. [51 the features do diminish as the number 
of planes of hit information increases and they do 
reverse when the field reverses. Further, there is a 
large dip at low momenta in the DF data set that 
is not seen in the DR data set which also becomes 
negligible as the number of planes of hit informa- 
tion increases. We therefore adopted the MIC cut 
at 60— planes to minimize these systematics. The 
cut was placed at 60 planes because the addition 
of more planes of information did not appreciably 
reduce the systematic error while it does reduce the 
statistics for the number of events passing the cut. 

(2b) "BdL" analysis - a track was required to have an 
integral field strength of BdL > 12 Tesla-meters 



(Tm), where BdL is a measure of the perpendicular 
magnetic field Bperp traversed by the track. For 
this analysis a variable BdL was defined as 

nend 

BdL = / Bperpir) dL, (1) 

J beg 

where Bperp — \ B{r) x Ti\ is the component of the 
magnetic field perpendicular to the track direction 
"n at a given point along the track path, r is the 
distance from the detector center axis, dL is the 
differential pathlength element along the track, beg 
is the point at which the muon enters the detector 
and end is the point at which it either exits the 
detector or stops in the detector. For our purposes, 
the track trajectory was approximated as a straight 
line running from the start to the end of the track. 
Since the track length of a long track is 43% in the 
steel, the cut at 12 Tm corresponds to a pt kick 
from the magnetic field of 1.5 GeV. 

The charge ratio as a function of the variable BdL 
is shown in Fig.[71 For low values of BdL, the mea- 
sured charge ratio approaches unity, as expected in 
the case of random charge determination. As the 
integrated magnetic field increases, the charge ra- 
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FIG. 7; The charge ratio A'^^+ /N^- as a function of BdL. 
The errors shown are statistical. Superposed is the fit to a 
constant charge ratio A'^+ l^ij.- for BdL > 12. 



tio increases to a plateau value where the charge 
misidentification is highly suppressed. The plateau 
value is reached at BdL = 12 Tm, as shown in 
Fig. [71 Above this value the charge ratio as function 
of BdL is consistent with being constant. 

The effect of these cuts on the DF and DR data samples 
is given in Table 1. 

As can be seen from Table HI tight cuts are required to 
minimize systematic errors in charge sign identification 
for the cosmic muon data sample. Since the Forward 
(DF) and Reverse (DR) data samples behave similarly 
with respect to the total suite of cuts, we can use both 
data sets in subsequent analyses, even though there are 



overall systematic differences in their iV„ 



-/Nfj_- ratios. 



IV. MUON CHARGE RATIO AT THE MINOS 
FAR DETECTOR 



In the sections below we determine the muon charge 
ratio, /N^- , as a function of the reconstructed muon 
momentum, p/it, by using the magnetic field of MINOS. 



A. Measurement of the Muon Charge Ratio 
Underground 

In Fig. [8] we show the charge ratio for events with re- 
constructed tracks that pass all cuts in Table |T] as a func- 
tion of p/it, the reconstructed momentum of the muon. 
In this and the figures that follow, p/^t is labeled "Mo- 
mentum" . In Fig. [8] the charge ratio is shown separately 
for (a) the DF and (b) the DR data sets; for these dis- 
tributions the MIC cut was used. Superposed onto the 
data are fits to a constant charge ratio. Fig. [5] clearly 
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Momentum (GeV/c) 




\) 50 100 150 200 250 

Momentum (GeM/c) 



FIG. 8: The charge ratio N^+fN^^^- as a function of fit mo- 
mentum for (a) the DF data set and (b) the DR data set. 
For this figure, the data set was selected using the MIC cut. 
Superposed on both are the fits to a constant charge ratio 



shows the systematic differences in the charge ratio mea- 
surements between Forward and Reverse field running. 

A method to cancel geometrical acceptance effects and 
alignment errors is discussed in [3, [3]. If Ai is the 
acceptance for /i+ and A2 is the acceptance for /i^ in 
the Forward field direction, then the acceptances in the 
Reverse field direction are Ai for fj,^ and A2 for . We 
can thus write two independent equations for the charge 
ratio in which the geometrical acceptances cancel: 



ra = {N^+/t)DF/{N^-/t) 



DR, 



and 



n = {N.+ /t)DB/iN.-/t) 



DF, 



(2) 



(3) 



By eliminating the Forward and Reverse live times, t^p 
and t]jji, between the two equations, we obtain a mea- 
surement of the mean charge ratio, r, in which both ge- 
ometrical acceptance and live time biases cancel, 

r=[rax ^ [(A^p+/A^p-)df x (iV^+/Af^- )z5fl]^/^ 

(4) 

Eq. ([U shows that it is the geometrical mean of the two 
independent charge ratio measurements that corrects for 
geometrical acceptance. 

Another class of systematic uncertainties, those that 
vary linearly with time, can also be cancelled by care- 
ful selection of the data analyzed. If there are systematic 
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TABLE I: Summary of the Cuts Applied 





DF" 


DR" 


-ff- c Veil La UfzIiUit; (^LlLo 


N=29.0 X 10'' 


N=8.9 X 10" 


cut 


Fraction Remaining 


No Cuts 


1.0 1 1.0 


Pre- Analysis Cuts: 




1. no reconstruction 


0.790 


0.832 


2. multiples 


0.733 


0.776 


3. coil status 


0.730 


0.772 


Analysis Cuts: 




1. 20 plane cut 


0.554 


0.585 


2. 2m track length cut 


0.551 


0.582 


3. fiducial cut 


0.534 


0.565 


4. fit quality cut: xj itter / ndf < 1.5 


0.427 


0.452 


Charge-sign quality cuts 




1. > 2.2 


0.141 


0.147 


2a. M/C cut 


0.048 


0.050 


2b. BdL cut 


0.033 


0.031 



" DP — cosmic data set, forward field 
' DR = cosmic data set, reverse field 



effects that depend on live time in a linear way, these sys- 
tematics can be cancelled by using Forward and Reverse 
data sets obtained during equal intervals of live time. In 
this analysis we make two independent measurements of 
the charge ratio using data sets constructed in this way. 
We first divide the Reverse data into two sets, DRi and 
DR2, each of which has live time equal to ^cis/S. We 
then pair DRi with DFi, a Forward data set also with 
live time tun/ 2 which falls at the end of the first period 
of Forward running (c.f.. Fig. S]), in Eq. These data 
result in one independent measurement of the charge ra- 
tio, ri, 

n = [{N^+/N^-)df^ X {N^+/N^-)drA^ (5) 

The second measurement comes from pairing DR2 with 
DF2, a Forward data set with live time 1 011/2 which 
comes at the beginning of the second period of Forward 
running (c.f.. Fig. d]), in Eq. These data result in 
second independent measurement of the charge ratio, r2, 

r2 = [{N^,^/N^-)df, X {N^^/N^-)dr,]^ (6) 

For these measurements, the data sets are all of length 
tDj^/2. This analysis relies on systematic errors domi- 
nating over statistical errors since the method limits the 
sample sizes. 

In sections IIV Bl and IIV CI the analysis just described 
is applied to the data sets generated with the MIC and 
BdL analyses. 

B. MIC Analysis 

In Fig.[5]we show the charge ratio /N^- for the two 
data sets ri and r2 as a function of p fa ■ The two inde- 
pendent measurements of the charge ratio are consistent 
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PIG. 9: The charge ratio N^+/N^- as a function of fit mo- 
mentum, Pfit- For this figure, the data set was generated 
using the MIC cut. The errors shown are statistical. Super- 
posed is the fit to a constant charge ratio A''^j+ /N^~ . 



with one another, which suggests we have significantly 
reduced the systematics seen in Fig. [5] We have used the 
data for ri and r2 independently in a fit to a constant 
charge ratio over the range < p^j^ < 150 (GeV/c). 
In this fit, each fit momentum bin was assumed to have 
two independent measurements of the charge ratio. The 
results of this fit give N^+/N^- = 1.370 ± 0.003 with 

/ndf = 1.15 for 29 degrees of freedom. The best fit to 
a constant charge ratio has been superposed onto the ri 
and r2 data in Fig. [Hland is labeled r. 

In Figs. [TUl and [TTl we show the charge ratio as a func- 
tion of zenith angle and azimuth. The analysis for both 
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TABLE 11: Number of Events in MIC and BdL Samples 
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FIG. 10: The charge ratio N^+/N^- as a function of cos^, 
where 6 is the zenith angle. The errors shown are statistical. 
Superposed is the fit to a constant charge ratio + /A - . 
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FIG. 11: The charge ratio A^+ /A^- as a function of az- 
imuthal angle (j). The errors shown are statistical. The gaps 
are due to acceptance effects resulting from the planar nature 
of the detector. Muons with azimuths < 60° have had 360° 
added to their azimuth so that there is a continuous distribu- 
tion of muons from the North (270° — 420°). Superposed is 
the fit to a constant charge ratio A^^+ /N^- . 



horizon coordinates follows the analysis for the charge ra- 
tio as a function of fit momentum. There is no evidence 
for an angular dependence of the measured charge ratio. 

From Figs. [9l [TOl and [Til we find small differences in 
the charge ratio depending on the parameter used to bin 
the data. Some of these differences could result from 
the particular binning chosen for the different parameter 
plots. To remove differences due to binning, we computed 
the charge ratios in Eq. (O and Eq. ^ directly for all 
data with p/^t < 250 GeV/c. Using the MIC selection 





MIC 


BdL 












DFl 


132,905 


94,792 


75,360 


52,897 


DRl 


128,789 


96,380 


66,581 


50,378 


DF2 


133,382 


94,434 


73,889 


51,126 


DR2 


125,526 


93,802 


67,520 


51,057 



in the charge ratio analysis, we find the bin free charge 
ratio is 



1.372 ±0.003, 



(7) 



where the error is statistical. The number of events used 
in the computation of this ratio are given in Table ITTl 



C. BdL Analysis 

Using the BdL selection in the charge ratio analysis, we 
find the bin free charge ratio for all data with p/^t < 250 
GeV/c is 



1.377 ±0.004, 



(8) 



where the error is statistical. The number of events used 
in the computation of this ratio are given in Table [Til 



D. Systematic Uncertainties 

Using the analysis technique described above many 
systematic uncertainties cancel. However, there are 
residual systematic uncertainties which do not cancel. 
These systematics can be separated into two classes: bias 
errors and random charge identification errors. Bias er- 
rors are those that lead to misidentifications of charge 
sign and cancel with a high degree of precision when 
combining forward and reverse field data. The extent to 
which these errors do not cancel is a measure of the mag- 
nitude of these systematics. Randomization errors lead 
to random misidentifications of charge that do not can- 
cel in Eq. (|3|). The magnitude of these latter systematic 
errors have been determined by Monte Carlo simulation. 
An example of a process that leads to randomization er- 
rors is multiple scattering, which can make a straight 
(high momentum) track appear curved to the track fitter. 
The charge misidentification of these tracks is indepen- 
dent of charge sign. A second source of randomization 
error comes from the inclusion of spurious hits in the 
track fit from bremsstrahlung, cross talk (signal appear- 
ing to come from a channel adjacent to the one hit), or 
the incorrect assignment of optically summed hits to a 
plane. Since these hits are likely to fall on either side 
of the track with equal probability, the incorrect charge 
sign assignment again occurs with equal probability. 

We examine three sources of systematic uncertainties. 
The first two are bias errors and the third is due to ran- 
domization errors. 
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1. Combining Forward and Reverse Data. Here we 
estimate the residual bias associated with focusing, 
detector acceptance, the magnetic field map, and 
misalignments. These cancel in principle when the 
forward and reverse data sets are combined with 
Eq. ((!]). As shown in Fig. [9l the method mostly 
suppresses the unphysical structures apparent in 
the charge ratio. We take (ra — rb)/2 = 0.009, 
or the residual difi'erences in the charge ratio due 
to incomplete cancellation of bias systematic errors 
remaining after the field has been reversed, as a 
measure of the systematic error associated with our 
procedure for combining forward and reverse data. 
The data sets used to calculate ra and rf, are the 
same as those used in the computation of ri . 

2. Sliding Window. This bias error determines how 
stable our determination of the charge ratio re- 
mains with respect to time. As described in IIV Al 
systematic errors which grow or decrease linearly 
in live time can be reduced by combining Forward 
and Reverse Field data sets in Eq. ^ which have 
equal live times IbrI^- These determine the charge 
ratios r\ and r2 . We have recomputed r\ and r-i by 
successively sliding a window of width IdbJ^ to al- 
ternative times in the first and second periods of 
Forward Field running. In this way we found a 
range of ±0.005 for the charge ratio and this range 
was adopted as the systematic error. 

3. Random Charge Identification. Random charge 
identification errors or 'randomization' errors are 
those that lead to events getting random charge 
sign assignments. Since the charge ratio iV^+ /N^- 
is greater than one, more positive than nega- 
tive muons will be misidentified and the measured 
charge ratio will decrease toward unity. Or equiva- 
lently, randomization errors always result in a mea- 
surement of the charge ratio that is lower than 
its true value. Consequently, randomization errors 
lead to a one-sided (positive) error on the measured 
charge ratio. 

We estimated the magnitude of the randomization 
systematic error with Monte Carlo simulation and 
the BdL analysis. We use the measurements of 
the charge ratio with successively tighter cuts on 
the BdL parameter (shown in Fig. [7]), a process 
that systematically reduces the randomization er- 
ror. Using the trend from Monte Carlo in a fit, 
we estimated the charge ratio without randomiza- 
tion errors. The difference between the fitted value 
of the charge ratio without randomization errors 
and the averaged (MIC/BdL) value of the charge 
ratio is 0.007. We adopt this value for the system- 
atic error due to random charge identification. 

Table IIIII summarizes these systematic uncertainties. 
The total systematic uncertainty was computed as the 
quadratic sum of individual uncertainties. 



TABLE III: Systematic Uncertainties in the Charge Ratio 



Source 


(^ST/s£ ) i 


1. Combining Forward/Reverse 

2. Sliding Window 

3. Randomization 


± 0.009 
± 0.005 
+0.007 




(+0.012, -0.010) 



E. Muon Charge Ratio Underground 

For the charge ratio underground as measured by MI- 
NOS, we take the average of the charge ratio obtained 
from the MIC and BdL analyses, and use the systematic 
uncertainties listed in Table Hill to give the muon charge 
ratio underground: 

r = 1.374 ± 0.004 (stat.)^!^:"^^ (sys.). (9) 



V. THE MUON CHARGE RATIO AT THE 
SURFACE 

To infer the muon momentum at the surface from the 
momentum measured underground requires knowledge 
of the rock overburden above the MINOS far detector. 
However, the rock overburden above the detector is of 
non-uniform composition with bands of iron formation 
embedded in Ely- Greenstone [l^. The topography of the 
surface above the MINOS far detector is also not level but 
rather has surface elevations that vary from 630 to 720 
m over the angular region of interest [17[ • Since the vari- 
ations in the composition of the rock overburden are not 
known directly, the technique used here is to normalize 
the data to a "world survey" of vertical muon intensity 
data [13. Thus, it is possible to derive a value of the 
slant depth for each solid angle bin. 



A. Projection back to the Surface 

1. Measured Vertical Muon Intensity at the MINOS Far 
Detector 

In this analysis we use events with good charge sign 
and momentum reconstruction, that is, those events that 
pass the cuts in Table|T]up to and including the M/C cut. 
As seen in TableHH this data set has the largest statistics. 
The selected events were first separated into bins of equal 
solid angle, AO = Acos(6') x A(f> ^ 0.02 x 0.10 sr. In 
each solid angle bin j, the vertical muon intensity was 
computed according to 



"'^ T {ejAjAn/cos0j)' 
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where T is the hve-time; Nj is the number of single muons 
in bin j; m is a multipUcative factor, assumed indepen- 
dent of slant depth and direction, that accounts for muon 
multiplicity; = e(cos 9 j,4'j) and Aj — A{cos9j,4)j) are 
the efficiency and the projected area of the detector, re- 
spectively, as a function of zenith and azimuthal angles; 
and the cos{9j) factor corrects for the muon intensity 
zenith angle dependence at the surface to a good ap- 
proximation in the momentum range relevant here. The 
intensity is converted to vertical intensity to facilitate 
comparison with the world survey data. Each measure- 
ment of is an independent measurement of the ver- 
tical muon intensity in direction [cos 9 

The efficiency of the far detector, ej = €{cos 9 for 
the reconstruction of single muon tracks was computed 
with Monte Carlo generated muons. For this calcula- 
tion, we generated a sample of 1.2 x 10^ cosmic muons 
by Monte Carlo simulation. Each event was generated 
by first choosing an arrival direction from the zenith an- 
gle dependence of the muon flux parameterization on the 
surface [l^ and then associating this direction with the 
overburden [ g/cm^] in the Soudan 2 slant depth map 
[l3]. The energy of the event was selected from the sur- 
face cosmic ray muon distribution [l^ . With the energy 
and overburden, the event was tested to see whether it 
penetrated to the detector [205 • The events that survived 
were placed on an imaginary box positioned around the 
detector [2lj and then propagated through the detector 
with the GEANT3 simulation of the MINOS detector. 

The standard method for computing the statistical er- 
ror in the efficiency is to consider the application of the 
cuts to be a binomial process. With the large number 
of events in our Monte Carlo sample, the statistical er- 
ror is much less than 1%. However, our Monte Carlo 
does not include multiples, a fair sample of demultiplex- 
ing failures, or events arriving when the magnetic field is 
off and these effects outweigh the statistical errors. Long 
experience with our MINOS Monte Carlo suggests that 
the uncertainty on our computation of e is of the order 
of a few percent. 

The projected area of the MINOS far detector in di- 
rection {cos9j, Aj — A{cos9j,(j)j), was computed 
by first finding the unit vector along this direction, 
— (cos 9j ,4'j)j then defining the normal for 
each of the ten surfaces of the MINOS far detector, fik, 
where fc = 1 — 10. The projected area is then given by 



vertical intensity to account for muon multiplicity, the 
correction factor m in Eq. (jlOp . Using a day of MINOS 
far detector data, we find that m — 1.04. Corrections 
were not made for the lateral distribution of multiple 
muons over the finite size of the MINOS detector. These 
corrections are quite labor-intensive to compute and we 
conservatively estimate them to be ^10%. The effect 
of these uncertainties on the analysis are considered in 

2. Determination of the MINOS Slant Depth Map for 
Standard Rock 

For each solid angle bin j in direction (cos 9j , (t>j ) we 
computed a value for the slant depth, Xj , or the column 
of rock from the MINOS cavern to the surface in units 
of m.w.e. (where 1 m.w.e. = 10^ g/cm^) by equating our 
measured vertical muon intensity to Crouch's all-world 
average vertical muon intensity |18j . This parameteriza- 
tion represents the integral of muons of all energies at 
the surface that can reach the detector through a rock 
depth Xj. In particular, we varied Xj until the Crouch 
parameterization and our measured vertical intensities 
agreed. Some solid angle bins at the edge of the accep- 
tance had too few events to reliably determine the slant 
depth; these solid angle bins and their events were re- 
moved from further analysis. One solid angle bin located 
near the zenith and on the edge of the detector accep- 
tance had just enough events for the slant depth of that 
bin to be calculated. However, the calculated slant depth 
was several hundred m.w.e. less than the slant depth of 
neighboring solid angle bins. Given the small acceptance 
of the bin and the peculiar slant depth derived for it, this 
solid angle bin and its events were also removed from the 
analysis. 

Since the Crouch parameterization is given in terms of 
"standard rock" , the MINOS slant depth map computed 
here is in terms of standard rock. In Fig. [12] we show 
the vertical muon intensity in the range 2000 m.w.e. < 
X < 5000 m.w.e. made using the MINOS slant depth 
map with the Crouch parameterization superposed. The 
two distributions coincide as they must. 

3. Vertical Muon Intensity for Soudan Rock 
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Aj = A{cos9j,(l}j) = J^^^t" ■ nk)5'fc, (11) 

k=l 

for all (n^ • fik) > and where Sk is the area of the 
fcth surface of the MINOS far detector. Our compu- 
tations show that the total acceptance for the MINOS 
far detector to single atmospheric muons is {eAft) — 
J2j(.ejAjAn) = 1.3 X 10^ cm2 sr. 

To compare our results to those of other underground 
experiments jl^, we make a correction to our measured 



Muons of energy E^^^o at the Earth's surface lose energy 
[Tgj as they traverse a slant depth X through the Soudan 
rock to the MINOS detector according to 

-^^a{E^) + b{E^)E^, (12) 

where the parameters a and b describe the energy lost by 
collisional and radiative processes, respectively. Eq. (fT2)) 
assumes continuous energy loss and does not account for 
fluctuations [l^]. The energy loss parameters for stan- 
dard rock, (asjbs), as a function of energy are given in 
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TABLE IV: Rock Parameters 
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FIG. 12: The vertical muon intensity in the MINOS far de- 
tector hall. Shown as solid circles are the vertical intensity 
data which have been normalized to the Crouch all-world av- 
erage vertical intensity for standard rock [l^. Overlaid onto 
the solid circles is the Crouch average. Shown as x's is the 
vertical intensity for the MINOS far detector hall corrected 
for Soudan rock. 



[19|. At the detector, the energy of the muons, i?^, can 
approximately be related to E^^ by 19] 

^^^,0 = (^p + a/6)e^^ - a/b. (13) 

To convert the MINOS slant depths for standard rock to 
Soudan rock, we equate the minimum energy required to 
reach slant depth Xg of standard rock to the minimum 
energy required to reach the equivalent slant depth Xm 
for Soudan rock [23|, |2^, l25|. If (aM,&M) describe the 
energy loss parameters for the Soudan rock at MINOS, 
then 

XM = T^ln[l + (^)(^)(e^=^=-l)]. (14) 



The values of the energy loss parameters a and b depend 
on the average composition of the rock [2^. For the 
collisional term, 

, Zm , , -^s , _i 

an ^as{—){—) , 

Am As 



and for the radiative term 



hM = bs{-^){-^ 



y1 y2 



where Z is the atomic number and A is the atomic 
mass. The parameters for standard rock and the Ely- 
Greenstone rock at the Soudan Underground Laboratory 
[HI are shown in Table ITVl 

The vertical intensity using the Soudan rock map is 
shown in Fig. [121 Since {Z\^IAm) is larger for Soudan 
rock than for standard rock, a muon will lose more en- 
ergy traversing an equivalent column of Soudan rock and 
therefore the vertical intensity at MINOS will have a 
steeper slope, as seen. In Fig. [13] we use the MINOS 
slant depth map to plot the charge ratio as a function of 
the slant depth. Fig. [13] shows that the charge ratio has 
little dependence on slant depth. 
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FIG. 13: The charge ratio A^^j+ /A^^- as a function of slant 
depth in Soudan rock. The errors shown are statistical. Su- 
perposed is the fit to a constant charge ratio /A^- . 



J^. Projection back to the Surface 

The muon energy underground is obtained from the 
reconstructed momentum. To project the reconstructed 
muon energies underground, -B^ = cpju, back to the sur- 
face, we used Eq. (|13p with Soudan rock parameters and 
the MINOS rock map. In these projections, we assumed 
that the Soudan energy loss parameters varied with en- 
ergy in a manner similar to the energy loss parameters 
for standard rock in [l^. In Fig. [T?] we show the result 
of this projection for three slant depth bins used to com- 
pute the vertical intensity in Fig. 1121 This figure shows 
the distributions of surface muon energies, £^p,o, for these 
three bins, as well as the median energy, (-E'^,o)med, for 
each distribution. The width of the slant depth bins is 
100 m.w.e. and this width is the dominant contributor to 
the width of the surface energy distributions shown. The 
expected increase of surface energy with increasing slant 
depth is clearly evident. 



The Energy Dependence of the Muon Charge 
Ratio at the Surface 



Using Eq. we projected the muons from the MIC 
analysis back to the surface. For each successfully re- 
constructed muon underground, we use its pfa value, its 
slant depth X, and the Soudan values of (aM, ^m) to ob- 
tain the surface energy -E^,o- In these projections, we 
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TABLE V: Charge Ratio at the Surface 
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FIG. 14: The distribution of muon energies projected back to 
the surface for three slant depth bins of width 100 m.w.e used 
in the computation of the vertical intensity in Fig. [121 These 
projections were made with Eq. (|13|l and Soudan rock param- 
eters. The median value of the muon energy (£'fi.o)med on the 
surface for these three bins is shown. 
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FIG. 15: The muon charge ratio A^^+ jN^- at the Earth's 
surface. The errors shown are statistical. 



assumed that the Soudan energy loss parameters were 
independent of charge sign and they varied with energy 
in a manner similar to the energy loss parameters for 
standard rock in p^ . 

Once projected back to the surface, we sorted the 
muons into bins of width 0.25 TeV and then computed 
the charge ratio. The results are shown in Fig. [151 where 
the charge ratio iV^+ jN is plotted as a function of sur- 
face muon energy i?^,o- In Table IVl we give the charge 
ratio data in each energy bin, as well as the weighted 
mean. 

We have performed 1-parameter and 2-parameter fits 
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to the data in Table |Vl over the surface muon energy 
range 1.0 TeV < S^,o < 7.0 TeV. A fit to a constant 
charge ratio gives 



1.371 ± 0.003 (stat. 



\ +0.012 
/-O.OlO 



(sys.), (15) 



with /ndf = 63.2/67. The data are consistent with a 
charge ratio that is independent of energy. 
The 2-parameter linear fit gives 

Nf,+ /N^- = (1.354 ± 0.007) + (0.85 ±0.33) x 10"^^^^^ 

(16) 

where -E^.o is in TeV and x^/i^df — 56.6/66. 

There are two contributions to the systematic error 
on the slope: uncertainties in the calculation of the en- 
ergy scale Efj^^Q and errors on the determination of the 
slope due to randomization. Errors in the energy scale 
are mostly due to uncertainties in the slant depth map. 
Using Monte Carlo methods to study these slant depth 
uncertainties, the errors are found to be ~10%. Cal- 
culations show that errors in the surface muon energies 
resulting from systematic uncertainties in the slant depth 
map of this order are 15-20% at 2100 m.w.e. and 25% at 
4000 m.w.e. In a second test, the Soudan 2 slant depth 
map fv^ was substituted for the MINOS slant depth map 
and the surface energies recomputed. The differences in 
the surface energies determined with these two maps are 
again approximately 20%. We therefore estimate the sys- 
tematic error on the energy scale to be ± 20%. The 
uncertainty in the energy scale does not affect the signif- 
icance of the slope. 
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Systematics due to randomization do affect the signifi- 
cance of the slope determination. For long track lengths, 
the tracks are better reconstructed, the charge misidenti- 
fication is smaller, and the measured charge ratio system- 
atically rises. Thus randomization can mimic and cover 
up a rising dependency of the muon charge ratio on muon 
surface energy. To estimate the size of this effect, we have 
separated the data into six subsets with BdL values from 
10 Tm to 15 Tm. As expected, we found these data sets 
to have different values of the charge ratio due to differ- 
ent but fixed amounts of randomization at the level of 
approximately 2%. In contrast, the values for the slope 
of the charge ratio versus surface energy were similar. 
This randomization-free value for the slope parameter 
was 30% smaller than the value from the two parameter 
fit. We take this difference to be the systematic error on 
the slope, 

slope ^ 0.85±0.33(stat)±0.26(syst) x 10"^ TeV"^ (17) 

Adding the errors in quadrature shows that the slope 
differs from zero by two sigma. An alternative fit using 
a linear dependence on log i?^ yields a similar result. 



VI. DISCUSSION 

In order to obtain the charge ratio results shown in 
Fig. 1151 each muon in our sample was projected back to 
the surface with Eq. (fT5)) using the Soudan energy loss 
parameters (omj^a/) for both /i+ and fi^ . A theoreti- 
cal complication to this procedure is the possibility that 
the energy loss parameters a and b are different for the 
two charges. For ionization losses, it was pointed out by 
Fermi [26?| that the electrons in matter would introduce 
small differences in energy loss. Calculations by Jackson 
and McCarthy [131 confirmed that negative particles lose 
energy at a slower rate, with the difference dropping from 
tens of percent at MeV energies to about 0.3% in the GeV 
range. These calculations were subsequeiitly verified ex- 
perimentally both at MeV energies [H, [23, [13] and in 
the GeV range [3l|. The approximations used in ^t\ . 
however, break down when going to even higher energies 
and so more exact numerical methods are needed. A cal- 
culation at TeV energies can be found in Jackson [32| 
that shows a small 0.15% increase in the ionization loss 
for /i"*" over For radiative energy loss, which involves 
the parameter 5, the difference between /i+ and /i^ is 
much smaller and falls with energy [33|. In the extrapo- 
lation used to obtain the results below, we assumed the 
same energy loss function for both charges. 

The projections of our data back to the surface, plot- 
ted in Fig. [15] as a function of surface energy, yield a 
charge ratio significantly higher at few TeV than those 
measured by others at surface energies below 300 GeV 
[1, [3. This rise in the charge ratio at TeV energies is, 
however, expected as the result of the increasing contri- 
bution of kaons to the cosmic ray muon flux at these 



energies and the greater likelihood for kaons to decay to 
^+ than for pions to decay to /x+ [20| . 

We use a qualitative model of the charge ratio to show 
that the rise in the charge ratio at TeV energies seen in 
Fig. [15] is consistent with this expectation. Gaisser [23| 
and Gaisser and Stanev {19] give an expression for the 
muon intensity at the surface as a function of the muon 
energy and zenith angle. It has contributions from both 
pion and kaon decay and it comes from folding the mea- 
sured spectrum of cosmic ray primaries with the kinemat- 
ics of pion and kaon decay. In our model we assume this 
expression holds independently for both /i+ and fi^ . In 
addition, this model assumes that charm production can 
be neglected and that pion and kaon interaction lengths 
are independent of charge at these energies. Let /jr+ be 
the fraction of all pion decays with a detected muon that 
is positive. Then (1 — /Tr+) is the fraction of all pion de- 
cays with a detected muon that is negative. Similarly, let 
fx:+ be the fraction of all kaon decays with a detected 
muon that is positive and (1 — fK+) be the fraction of all 
kaon decays with a detected muon that is negative. The 
muon charge ratio /N^- can then be written 



N„ 



N„ 



I.IE,, 



+ 



1.1E__-L cos 6 



(1 



115 GcV 
/.+ ) 



1 



850 GcV 

0.054(1 



1 + 



I.IE _ 



115 GoV 



1 



I.IB _ • 



850 GcV 



(18) 



The simplest assumption to make in this model is that 
/^+ and fx+ are independent of energy. Although this 
assumption neglects many physical processes that may 
play a role at these energies, it is an assumption that is a 
reasonable choice to qualitatively describe our results: a 
rise in the charge ratio from a plateau at a few hundred 
MeV to a second higher plateau at a few TeV. This rise 
has already been seen in the results from the CORT cos- 
mic ray Monte Carlo [H and the models of Lipari [sl ]. 
Below we test this simple model with the MINOS data. 

We used Eq. (Tip with the MINOS data set and the 
L3-I-C data set ;3| to find the values of /7r+ and fx+ 
that best describe these two data sets. We used these 
data because they have the angular information needed 
for the fit. We found the best fit values for and /^+ 
with a grid search over {fTr+,fK+) parameter space. At 
each point in the space a statistic compared the mean 
charge ratio weighted by solid angle with the model pre- 
dictions represented by Eq. ((TS]). The charge ratio values 
in each bin had uncertainties given by the quadratic sum 
of the statistical and systematic errors. The min- 
imimum is found at /^+ = 0.55 and = 0.67, with 

/ndf ~ 1. In Fig. we have superposed this 'ttJC' 
model onto the MINOS and L3-fC data sets [1]. The 
qualitative results of the model are that the observed 
rise in the muon charge ratio can be explained by the in- 
creasing importance of kaon decays to the muon charge 
ratio as the energy increases from 0.3-1 TeV and that 
values of and fj^+ that are independent of energy 
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FIG. 16: The nK model discussed in the text superposed onto 
the MINOS and L3+C data sets. The MINOS charge ratio 
data are from Table FVl with a systematic error of ±0.011 that 
has been added in quadrature to the statistical error for each 
data point. The L3+C data are taken from [^. 
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FIG. 17; A compilation of muon charge ratio data from 0.1 
to 7 TeV. The MINOS data have been taken from TableEl 
Other data: L3+C Baxendale 35], CosmoALEPH [sj 
Matsuno [3, and Rastin 37]. The ttK model is superposed. 



are sufficient to describe the MINOS and L3+C data. 

Fig. 1 171 shows the ttK model superposed onto a compi- 
lation of data from the literature, as well as the MINOS 
and L3+C data. The additional data also support the 



results of this simple model. 



VII. SUMMARY 

The analysis presented here can be separated into three 
parts. First, we computed the muon charge ratio un- 
derground. To minimize residual systematic errors, we 
combined data with both foward and reversed magnetic 
field running configurations. When combined, two inde- 
pendent analyses (MIC/BdL) give a muon charge ratio 
as measured underground of: 



1.374 ± 0.004 (stat.)+|^:^io (sys.). 



Second, using a map of the Soudan rock overburden, the 
muon momenta were extrapolated to their corresponding 
values at the surface, spanning the energy range from 1 
to 7 TeV. Within this range of energies at the surface, 
the MINOS data are consistent with the charge ratio be- 
ing energy independent at the two standard deviation 
level. The charge ratio as measured by MINOS is signif- 
icantly higher than measurements by other experiments 
at surface energies below 300 GeV. Finally, we used MI- 
NOS and L3-I-C data in a simple model that attributes 
the rise in the charge ratio with energy to the increas- 
ing contribution of kaon decays. Fitting the data to the 
model gives results that are consistent with this picture. 
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